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ABSTRACT 

Far- infrared spectra of laboratory silicates are normally obtained at room temperature 
even though the grains responsible for astronomical silicate emission bands seen at 
wavelengths > 20 /xm are likely to be at temperatures below ~ 150 K. In order to 
investigate the effect of temperature on silicate spectra, we have obtained absorption 
spectra of powdered forsterite and olivine, along with two orthoenstatites and diopside 
clinopyroxene, at 3.5±0.5 K and at room temperature (295±2 K). To determine the 
changes in the spectra the resolution must be increased from '--^ 1 to 0.25 cm~^ at both 
temperatures since a reduction in temperature reduces the phonon density, thereby 
reducing the width of the infrared peaks. Several bands observed at 295 K split at 
3.5 K. At 3.5 K the widths of isolated single bands in olivine, enstatites and diopside 
are ^ 90% of their 295 K-widths. However, in forsterite the 3.5-K-widths of the 31-, 
49- and 69-/im bands are, respectively, 90%, 45% and 31% of their 295 K widths. Due 
to an increase in phonon energy as the lattice contracts, 3.5-K-singlet peaks occur 
at shorter wavelengths than do the corresponding 295-K peaks; the magnitude of the 
wavelength shift increases from — 0.2 /j,m at 25 /im to ~ 0.9 /xm at 80 /im. In 
olivines and enstatites the wavelength shifts can be approximated by polynomials of 
the form ax + hx^ where x = Xpk (295 K) and the coefficients a and b differ between 
minerals; for diopside this formula gives a lower limit to the shift. Changes in the 
relative absorbances of spectral peaks are also observed. The temperature dependence 
of Xpk and bandwidth shows promise as a means to deduce characteristic temperatures 
of mineralogically distinct grain populations. In addition, the observed changes in band 
strength with temperature will affect estimates of grain masses and relative mineral 
abundances inferred using room-temperature laboratory data. Spectral measurements 
of a variety of minerals at a range of temperatures are required to fully quantify these 
effects. 
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1 INTRODUCTION 

Mid and far-infrared spectra obtained with the Infrared 
Space Observatory Short- and Long- Wavelength Spectrom- 
eters (SWS and LWS, respectively) have revealed emission 
bands which have been associated with crystalline silicate 
dust. The environments include comets, young stellar ob- 
jects and oxygen-rich dust in outflows and disks associated 
with late-type stars and planetary nebulae (e.g. Crovisier 
et al. 2000, Malfait et al. 1999, Waters et al. 1996, Sylvester 
et al. 1999 and Cohen et al. 1999, respectively). Simple spec- 
tral fits indicate that much of the optically thin silicate 
emission occurs at temperatures in the 50-100 K range (e.g. 
Sylvester et al. 1999). However, previous laboratory stud- 



ies (e.g. Mennella et al. 1998, Henning & Mutschke 1997, 
Day 1976, Agladze et al. 1996) have shown that the opti- 
cal properties of silicates vary between 295 K and 10 K. In 
MgO the frequency and FWHM(i/) of the 24.5 fim (295K) 
TO mode (Jasperse et al. 1964) vary nearly linearly with 
temperature between 7.5 and 1950 K (Kachare et al. 1972). 
A similar study of the effect of 423-873K-temperatures on 
3-12^m PAH bands by Joblin et al. (1995) revealed changes 
in band structure, bandwidth and frequency which indicate 
that free (as opposed to condensed) PAH molecules carry 
the astronomical 3.3/im band. 

Due to a lack of laboratory data, both simple black body 
fits and radiative transfer models of these environments have 
used silicate laboratory spectra obtained at room tempera- 
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ture. In addition to fitting ISO data, spectra obtained for a 
range of samples and temperatures will also be required for 
comparison with data from SIRTF, FIRST and the NGST. 
To begin to quantify the effect of temperature on silicates 
we compare spectra of two olivines, two orthoenstatites and 
diopside clinopyroxene at 295 K and 3.5 K. 



1.1 Choice of Spectral Resolution 

At cryogenic temperatures the vibrational states of a solid 
are substantially depopulated in comparison to their values 
at 295 K in accord with statistical thermodynamics. From 
the damped harmonic oscillator model, the width of a peak 
i s related to the lifetime betwee n phonon scattering events 
( IWooten 197^ [Hofmeister 1999| ). As the number of states 
decreases, the possibilities of scattering decrease, and the 
peaks must therefore sharpen. The vibrational bands ob- 
served in_cri;2S6nj£.J,0_^^J^00_jin^ crystalline sili- 
cates ( Day 1976 ; Henning fc Mutschke 1997 ) are sharper and 
deeper than those obtained at room temperature. Medium 
(~l-2 cm~^) resolution spectra (Henning & Mutschke 1997, 
Mutschke pers. commun. and Mennella et al. 1998) indicate 
that the bands of crystalline silicates shift in wavelength 
as the temperature is reduced, and that the magnitude of 
the shift increases with wavelength. If the spectral resolu- 
tion is insufficient, very narrow peaks are under- sampled, 
leading to rounded profiles and incorrect measurements of 
the wavelength shifts. Therefore, higher-resolution measure- 
ments are crucial in the far-IR where lower-resolution data 
(Mennella et al. 1998) showed that the isolated ^ 70 /xm 
(143 cm~^) peak of forsterite has a 24 K width of ~4 cm~^ 
and a ^2 cm~^ shift between 295 K and 24 K. We therefore 
present 0.25 cm~^ -resolution transmission spectra of crys- 
talline minerals from the olivine and pyroxene groups at 
3.5 K and 295 K. 
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Figure 1. Schematic three-dimensional structure of crystalline 
olivines. The lower diagram is a slice through the upper diagram 
along the h axis. Tetrahedra represent SiO^~ groups whose apices 
point out of the page (solid) and into the page (dotted). Ml and 
M2 represent the sites of divalent metal cations. 



2 LATTICE STRUCTURES AND 295 K BAND 
ASSIGNMENTS 

2.1 Olivines 

The general formula for the olivine solid solution series is 
(Mgi,Fei_a;)2Si04; the Mg^^ end-member is forsterite, the 
Fe^"*" end member is fayalite. Olivine minerals have the 
structure shown in Figure 0. In the figure each tetrahe- 
dron represents the Si04~ anion which consists of an Si 
atom surrounded by oxygen atoms at each corner of the 
tetrahedron; Ml and M2 represent the sites of the metal 
cations. In end-member forsterite both sites are occupied 
by Mg^+. In olivine some of the Mg^+ in sites Ml and 
M2 is replaced by Fe^"^; in this case there is no preference 
for substitution between the sites since these two ions have 
similar radii. The active modes in the 15-85-/.tm region are 
(after Hofmeister 1997): bending of the tetrahedral O-Si- 
O (A(295K) <, 20 ^m), rotation of the SiO*" tetrahedra 
(four bands between 21 and 26 /^m), translation of the Si04~ 
tetrahedra (~ 50 ^tm) and translation of the Mg'^"'' and Fe^"*" 
cations (several of the bands in the 21-30-/im region); two 
bands are associated with combined translations of the diva- 
lent cations and SiO*" tetrahedra (at ~ 36.5 and ~ 70 ^ra 
in forsterite; at ~ 94 and ^ 109 /xm in fayalite, Fe2Si04). 



2.2 Ortho- and Clino-pyroxenes 

Minerals in the pyroxene solid solution 
series {yigx ,Vey ,C&z)2'sihOfi where x + y -\- z = 1, belong to 
one of two crystal systems: orthorhombic, in which the three 
crystal axes are orthogonal, and monoclinic, which has two 
orthogonal and one inclined axis. The structure of orthopy- 
roxene is shown in Figure |^; silicate tetrahedra are joined 
together by shared oxygen atoms to form long chains along 
the c axis. Once again, Mg^""" and Fe'^'*' can replace each 
other in the Ml and M2 sites. The precise band assignments 
for the orthopyroxenes are unknown. However, comparison 
with the bands of the clinopyroxenes and olivines suggests 
that the bands longwards of 20 /im may involve translations 
of the cations, whereas bands at wavelengths shorter than 
20 /im involve the deformation of the tetrahedral and inter- 
tetrahedral 0-Si-O bonds. 

In clinopyroxenes similar to diopside (there are other 
types of clinopyroxene), the larger Ca'^"'" ion (r ~ 1.12A) 
occurs only in the M2 sites, whilst the smaller Mg^^ (r ~ 
O.72A) and Fe^+ (r ~ 0.78A) ions occur in both these and 
the smaller Ml sites. This changes the alignment of the 
tetrahedra, resulting in a crystal shape with an inclined a 
axis. Band assignments in diopsides are currently derived by 
comparison between force-constant calculations and labora- 
tory data. The active modes in the 15-85-^m region are mod- 
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Figure 2. Schematic three-dimensional structure of ortho- 
enstatite. The lower diagram is a slice through the upper dia- 
gram along the b axis. The b — c plane is characterized by chains 
of silicate tetrahedra joined by shared oxygen atoms. These are 
surrounded by metal ions (for clarity these are omitted from the 
lower pairs of tetrahedra in each chain). In the a — b plane Ml 
and M2 sites occur, respectively, at the apices and bases of the 
silicate tetrahedra. M2 sites are larger than Ml sites; therefore 
in diopside-like clinopyroxenes relatively large Ca^ 



ions can oc- 



cupy M2 sites but only smaller Mg^+ and Fe^+ ions may occupy 
the Ml sites. 



elled as (Tomisaka & lishi 1980) : bending of the tetrahedral 
and inter-tetrahedral 0-Si-O bonds shortwards of 30 /im 
at 295 K (at 23.2 and 29.6 /im these may be blended with 
Mg^''" translations) ; Mg^^ translations blended with transla- 
tions of the O-Si-O chain and rotations of the tetrahedra at 
32.3 /xm; Mg'^"'" translations blended with tetrahedral rota- 
tions at 34.3 /xm; Ca'^"'" translations blended with tetrahedral 
rotations at 40.4 fim; Ca^^ translations at 44.9 fj,m; finally 
a broad band at ~ 66 nm is produced by a zig-zag motion of 
bridging oxygen atoms in the silicate chain. However, force 
constant models are not backed-up by laboratory measure- 
ments, and the assignments are tentative. Measurements of 
the effect of isotopic or chemical substitutions are needed 
for reliable band assignments. 



3 THE EFFECT OF TEMPERATURE 



3.1 Changes in bond length 

If the bond length decreases, the energy of the associated 
phonon increases. (This physics is equivalent to the energy 
levels of a particle in a box; the effect is frequently observed 
in pressure studies, e.g. Ferraro 1984). Because frequency 
increases as pressure increases (or temperature decreases), 
the wavelength of the infrared peak must decrease. Tetra- 
hedral Si-0 bond lengths and bond angles are relatively 
incompressible in comparison to the spacings between the 
oxygen and metal ions and inter-tetrahedral oxygen atoms 
so there is little variation in the wavelength of the 10-/im 
stretching and 20-fim bending features with temperature 
(e.g. Day 1976). In contrast, the volumes of the divalent- 
cation sites are more compressible. Therefore one would ex- 
pect to see a greater change in the wavelengths of peaks 
associated with the cation sites. 



3.2 Changes in band occupancy 

The lifetime, F, of a phonon (the time between scatter- 
ing events) is given by F~^ oc 27rAi^(e2), where Az^(e2) 
is the full-width-half-maximum of the peak in the imagi- 
nary part of the dielectric function, £2, in frequency units. 
(Wooten 1972). When the bands are narrow and weak, the 
frequency for the atoms vibrating parallel to the propaga- 
tion direction of the radiation (the longitudinal optic; LO) 
is similar to the frequency of the atoms vibrating perpen- 
dicular to the propagation direction (the transverse optic; 
TO), thereby creating a symmetrical absorption peak whose 
shape approaches the Lorentzian shape of its counterpart in 
€2 (Hofmeister and Mao 2001). Absorbance full- width-half- 
maxima (FWHM) are related to dielectric widths; therefore 
the widths of the absorbance peaks should be roughly pro- 
portional to the dielectric widths in wavenumbers. 

At temperature, T, the occupation of phonon energy 
levels relative to the ground state is given by the partition 
function q = X]^o ^■''^P where k is the Boltzmann 

constant and it is assumed that the phonon frequency, v, 
is no t a function of temperature (this is not true; see Sec- 
tion 3.1). As the temperature is decreased the number of 



vibrational states decreases and the numbers of phonon- 
phonon scattering events decrease so that phonon scattering 
off the crystal sides and defects becomes more important. 
For phonon-phonon scattering, the lifetime is inversely pro- 
portional to the number of phonons available, and hence as 
the sample is cooled, states are depopulated, and the lifetime 
increases and peak width decreases. In practise, FWHMs 
are not often measured and the FWHM temperature de- 
pendence is not well understood (see Hofmeister, in prep). 



4 THE EXPERIMENT 
4.1 Sample Preparation 

Bulk samples of the minerals listed in Table |l| were bro- 
ken into smaller pieces with a hammer and then ground to 
powders by hand with a ceramic pestle and mortar. Grain 
size was measured by visual inspection in a binocular micro- 
scope at X 63 magnification. Small quantities of powder were 
mixed with petroleum jelly on 0.8-mm-thick polyethylene 
substrates. Preliminary room temperature measurements at 
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Table 1. Origin and composition of crystalline samples 



Group Mineral Composition 

Olivine Synthetic Forsteritt]^ Mg2Si04 

Olivine (var. peridot) (Mg1.77Feo.20 Alo.oiNio.oi)Si04 



Orthopyroxene Enstatite 

Bramble Enstatite 

Clinopyroxene Diopside 



(Mgi.8Feo.l6Cao.o2)(Sii.88Alo.l2)06 

{Mgi.7Feo.28Cao.oi)Si206 



Source 

Geophysical Laboratory, 
Washington, D. C., USA. 
Navajo Indian Reservation 
San Carlos, Arizona, USA 

Fiskenaesset, Greenland 
Bramble, Norway 



{Cao.992Nao.o28Mgo.94FeQ+^7Fe|^+j3Cro.oo2Alo.ai5)Si206 Dekalb, NY, USA 



Smithsonian sample R18682 



Synthesized from oxide powders by H. K. Mao 



Washington University indicated that the grain sizes and 
the numbers of grains on each substrate were sufficient to 
resolve most infrared bands. At Queen Mary, University of 
London (QMUL) the samples and substrates were cut down 
and placed in brass holders. The 295 K and 3.5 K measure- 
ments for each mineral were performed with a single sample 
and the same spectrometer. The only difference between the 
measurements was the positioning of the sample in the in- 
strument beam and the temperature of the sample. Hence, 
the substrate, the coverage of the substrate and the column 
density were identical for each measurement of a specimen. 



olivine and for thinner films of olivine, and the pyroxenes. In 
the thick films, bands shortwards of 45 //m were frequently 
saturated (these bands were excluded from the figures and 
the analysis). The olivine and enstatite thin films were too 
thin for bands longwards of 45 pim to be resolved. The Bram- 
ble enstatite and diopside films were of sufficient thickness 
for the majority of bands to be seen at both temperatures. 
Since the mass and grain size distribution of the silicate pow- 
der is unmeasured, the derived absorbances do not indicate 
mass density and the absorbance per unit area per unit mass 
cannot be estimated from these data at present. 



4.2 Cryogenics and Spectroscopy 

Spectra were obtained at QMUL using a Martin-Puplett po- 
larising Fourier Transform Spectrometer with a mercury arc 
lamp source and a liquid-helium-cooled bolometer pumped 
to 1.5 K. For the room temperature (295 K) measurements, 
the sample was placed in the output beam of the spectrom- 
eter; the 3.5-K tests were performed by putting the samples 
in a filter wheel in front of the detector within the 1.5-K 
bath in the cryostat. Since the sample was heated by in- 
cident radiation from the arc lamp, it was slightly warmer 
than the bath, at 3.5 ± 0.5K. The spectral resolution was 
selected by inspection of the interferogram; to resolve the 
fine structure of end-member forsterite at 3.5 K the resolu- 
tion must be increased to 0.25 cm"'^, whilst a resolution of 
2 cm~^ is adequate at 295 K. For consistency, all spectra at 
all temperatures were obtained at the higher resolution of 
0.25 cm"\ 

For each temperature, reference spectra of the open 
aperture and petroleum jelly on a polyethylene slide were 
obtained and then divided into the spectrum of the sam- 
ple to remove instrumental and substrate effects. For the 
3.5 K transmittance spectra, there is an uncertainty in the 
flux calibration. Therefore the 3.5 K data were normalised 
to the level of the 295 K transmittance spectra at 500 pm 
(20 cm~^) where the temperature dependence is small. The 
transmittance scale of the room temperature spectra is bet- 
ter than 1% and the uncertainty in the normalisation of the 
3.5 K data is ~ ±10%. The data were then converted to ab- 
sorbance units [A — — In (Transmittance)). Ripples in the 
spectra at wavelengths ^ 50 are an artifact of imper- 
fections in the thickness and the orientation of the samples 
and references relative to the beam. 

Data were obtained for thick films of forsterite and 



4.3 Astronomical Data 

For comparison of our laboratory data with SWS and LWS 
spectra from ISO, we present observations of CPD -56°8032, 
a planetary nebula whose ISO spectrum was published by 
Cohen et al. (1999). We have reduced SWS observation 
13602083 and LWS observation 08401538 and subtracted 
a polynomial continuum fit (with no physical significance) 
to show up the fine structure. Cohen et al. obtained a 
qualitative match to their spectrum with room-temperature 
forsterite (Mg2Si04), orthopyroxene and clinopyroxene, us- 
ing laboratory data by Koike et al. (1993) multiplied by 65K 
(forsterite, clinopyroxene) and 90K (orthopyroxene) black- 
body continua. Features at 43/im and 62^m were matched 
with crystalline water ice. The wavelengths and FWHMs of 
the silicate peaks in their fits are marked in Figures ^ and ^ 
by horizontal bars and vertical ticks. (It is beyond the scope 
of this paper to deduce the mineralogy of this source using 
our new laboratory data.) 



5 RESULTS 

5.1 Overview of the laboratory dataset 

The 15-85-/im absorbance spectra of the olivines (forsterite 
and olivine), orthopyroxenes (enstatite and Bramble en- 
statite) and a clinopyroxene (diopside) are presented in Fig- 
ure ^. Since it may be possible to estimate the grain size 
distribution and column density by electron microscopy at 
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Figure 3. Overview of data: 15 to 85 /im absorbance spectra of (a) olivines and (b) ortliopyroxenes (enstatite and Bramble enstatite) 
and a clinopyroxene (diopside) at 295 K (dotted) and 3.5 K (solid). Each 3.5 K spectrum is displaced by 0.5 units above the 295 K 
spectrum of the same sample. For thick films, bands shortwards of 45 fira were saturated; for thin films, bands longwards of 45 fira were 
too weak to be resolved; the Bramble enstatite and diopside films were of intermediate thickness so that most bands were resolved at 
both temperatures. 
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a later date, we present the measured values of absorbance 
for these samples^. 

In general, (i) the room temperature spectra look like 
low-resolution versions of their 3.5 K counterparts, even 
though the same spectral resolution was used for all these 
measurements; (ii) the relative absorbance of neighbouring 
peaks changes as the temperature is reduced; (iii) the wave- 
length of a peak decreases with temperature and the magni- 
tude of the shift increases with the wavelength of the band. 
The effect of temperature is most noticeable in the forsterite 
spectra, with the bands becoming particularly narrow and 
strong: this difference with degree of ionic substitution is 
real since this is the only mineral devoid of Fe^^ or other 
cations substituting for Mg^""" in the crystal structure. 



5.2 Effect of sample thickness 

Comparison of the spectra of the "thick" and "thin" olivine 
films reveals that thick samples (i.e. larger amounts of pow- 
der embedded in the petroleum jelly) are required to see the 
comparatively weak 50 and 70 ^m bands and this has re- 
sulted in the saturation of the most prominent features of 
forsterite and thick olivine at wavelengths < 40 /xm. Hence, 
saturated features have been excluded from the presented 
datasets. The room temperature spectra of the thin enstatite 
and thin olivine films contain no distinguishable features 
longwards of 45 /im. The Bramble enstatite and diopside 
films are of intermediate thickness since bands are not satu- 
rated and the far-infrared peaks are clearly distinguishable 
from the underlying ripple. 



5.3 Band splitting, band enhancements and 
changes in bandwidth 

The detailed structure in the 18-38-Atm and 38-80-//m re- 
gions are shown in Figures ^ and ^ respectively. Band 
splittings are seen shortwards of 38 /xm; whilst wavelength 
shifts, band enhancements and changes in the full-width- 
half-maxima (FWHM) of peaks are seen throughout the 18- 
80-/im range. 



5.3.1 Band splitting 

In both olivines the degenerate ~ 21 /im Mg^+ translation 
and Si04~ rotation splits into two bands when cooled to 
3.5 K (Figure |^. The ~ 34- /xm peak in forsterite is sat- 
urated, however this Mg^+ translation band is split in the 
3.5-K olivine spectrum and the 36.45-/xm 295 K Mg'^"'" trans- 
lation peak in forsterite is also split at 3.5 K. The band split- 
tings in forsterite are consistent with the peaks produced 
f rom different pol arisations occurring at the same frequency 
(Hofmeister 1997); at 295 K this is an unresolved doublet. 



Table 2. Peak wavelengths (Ap^ and FWHMs (AA) of selected 
bands. 



In the orthopyroxenes the band splitting is seen only in the 
three bands shortward of 23 ^m. Since the band assignments 
are undetermined in the enstatites, the cause of the splitting 



* This normalisation differs from that presented previously: 
Bowey et al. (2000) re-normalised pairs of 3.5 K and 295 K ab- 
sorbance spectra to the absorbance of the highest room temper- 
ature peak. 



Sample 



Forsterite 



Olivine 

Bramble 
Enstatite 
Diopside 



Apfc (295 K) 

^m 
31.27 ±0.02 
49.85 ±0.06 
69.74 ±0.12 
50.67 ±0.06 
73.20 ±0.13 
44.48 ±0.05 



AA (295 K) 

/xm 
0.29 ±0.02 
0.56 ±0.08 
1.3 ±0.1 

1.8 ±0.1 
2.76 ±0.05 



Apfc (3.5 K) AA (3.5 K) 



31.20 ±0.02 
49.35 ±0.06 
68.87 ±0.12 

50.21 ±0.06 
72.41 ±0.13 
43.94 ±0.05 



^m 
0.26 ±0.02 
0.25 ±0.08 
0.4 ±0.1 
1.4 ±0.1 
1.6 ±0.1 
2.49 ± 0.05 



40.37 ±0.04 2.59 ±0.05 40.19 ±0.04 2.24 ± 0.05 
44.93 ±0.05 1.77 ±0.05 44.62 ± 0.05 1.63 ± 0.05 
66.58 ±0.11 8.1 ±0.2 65.90 ±0.11 6.9 ±0.2 
The uncertainties in AA include the effects of spectral resolution, 
fringing and base-line estimates. 

is unknown. In the diopside clinopyroxene only the degen- 
erate Si-O-Si deformation at 20.91 /xm is split. 



5.3.2 Band enhancements 

The bands which are most strongly enhanced at low temper- 
ature, but which occur in both the room temperature and 
the 3.5 K spectra, are the unsaturated forsterite peaks, the 
25.44-(295 K) to 26.17-/im (3.5 K) olivine band, the 22.12- 
(295 K) to 23.00-^m (3.5 K) band of Bramble enstatite and 
the 25.22- (295 K) to 24.80-^im (3.5 K) band of diopside. 

Some peaks which do not result from obvious splitting 
(or changes in the substrate) are seen only in the 3.5-K spec- 
tra. Such features occur in the 3.5-K thin enstatite spectrum 
at 26.41, 29.17 and 34.06 ^m. Since similar bands are seen 
in Bramble enstatite spectra obtained at both temperatures, 
we suggest that these features are enhanced at 3.5 K. Other 
peaks are seen in Bramble enstatite but not in enstatite (e.g. 
the 295-K 22.12-, 31.98-, and 31.08-Aim bands); this is prob- 
ably due to the greater thickness of the Bramble sample and 
to differences in the metal content in the two minerals. Extra 
bands are seen in the 3.5 K spectrum of Bramble enstatite 
at 53.91 /im, and in diopside at 57.16 ^m. 



5.3.3 Changes in bandwidth 

Estimates of the full-width-half-maxima of peaks which are 
well-separated from adjacent peaks and well-resolved are 
given in Table ^; peaks which are unresolved doublets at 
295 K or are blended with stronger peaks have not been mea- 
sured. The widths of isolated singlet peaks provide insight 
into the population of the vibrational modes. FWHM val- 
ues were extracted with the IDL GAUSSFIT routine which ap- 
proximates the continuum local to the peak with a quadratic 
equatio n an d fits a gaussian to the peak. As described in 
Section 3.2, singlet peaks in the 3.5 K spectra are narrower 
than those measured at 295 K. With the exception of the 
forsterite peaks, bandwidths in the 3.5 K spectra are ~ 90% 
of their 295 K-widths, irrespective of the wavelength of the 
peak. The 31-, 49- and 69-/im bands of 3.5 K forsterite are, 
respectively, 90%, 45% and 31% of their 295 K widths; these 
bands arise from phonons along one direction within the lat- 
tice. 
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Figure 4. Detailed structure in the 18— 38-/xm wavelength range for the 295 K (dotted) and 3.5 K (solid) absorbance measurements, with 
the ISO spectrum of CPD -56°8032 (thick solid line) plotted for comparison. ISO observation: bold crosses indicate limits of spectral 
segments; horizontal bars with ticks indicate the FWHM and Apj. (respectively) of the features in the Cohen et al. (1999) model obtained 
using room temperature data: "F"=forsterite, "Op"=orthopyroxene, "Cp" =clinopyroxene) . Laboratory spectra: Lines and pairs of 
lines linking pairs of spectra respectively denote band shifts and band splitting. The styles of these lines indicate band assignments: 
solid— translation of Mg^+ and/or Fe'^^ cations, dashed- deformation of Si-O-Si bands, dot-das/i-rotation of SiO^" tetrahedra, dash- 
dot-dot-dot— translation of the silicate chain, dotted— the assignments of these bands are unknown; arrows denote bands seen only at 
3.5 K. Offsets in the y-axis are arbitrary; pairs of spectra are normalised to the same level; hence the relative changes in band strength 
are real for individual samples. 
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Figure 5. Band enhancements in tlie 38-80-/im wavelength range compared with ISO data for CPD -56°8032. The solid lines linking 
the laboratory curves indicate band shifts due to ionic translations, T{ion), except Z(-O-) which is a zig-zag motion of the bridging 
oxygen atoms between silicate tetrahedra; dot-dash lines indicate rotations of SiO^" tetrahedra, band assignments for dotted lines are 
unknown. Arrows indicate extra bands observed at 3.5 K. See Figure ^ for key to remaining symbols. 



© ??? RAS, MNRAS 000, §-0 



Temperature effects on the 15-85-^m Spectra of Olivines and Pyroxenes 9 



i 1 i i 1 i 1 i 1 i 1 1 i i i i 1 1 1 1 1 1 1 1 1 1 
: Forsterite (Thick) 


1 1 1 i 1 i i i i 1 i 1 i 1 1 i i i i 1 i i i i 1 1 F 




-'' *T(Si(U 
• T(Ug**) 
SiOli rotation 
t. Sl-O-Sl bend 
n^ig~zag of bridging 

' 


• 

1 1 . . . , 



i 1 i i 1 i 1 i 1 i 1 1 i i i 


1 i i i i 1 i i i 1 1 1 i 1 i 1 i i i i 1 i 1 i 1 1 i i i i 1 i i i i 1 i 1 i 1 

; ^\ 


". 1 . . , . 1 . , . , 1 . . . 


1 .... 1 .... 1 , ... 1 .... 1 1 .... 1 .... 1 ... ," 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




' 


.1 


Enstatite (Thin) 

' 



1 1 i i 1 1 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 i i 1 1 1 i 1 1 1 1 1 1 1 1 1 1 1 i i 1 i 1 i 1 


!^ i 






Bramble Enstatite 


' 








t \ 





I- 



Diopside 



^ Ml sites: Mg^ 



M2 sites: Ca^* 



I 20 25 30 35 HO U5 50 55 60 65 70 75 
Wavelength of Peak at 295K/Mm 

Figure 6. Shifts in the wavelengths of the peaks (A(295K) — 
A(3.5K)) versus the wavelength of the peaks at 295 K {A(295K)). 
Error bars are plotted where they are larger than the plotting 
symbols; empty circles indicate unknown band assignments; dot- 
ted lines indicate degenerate 295 K bands which split at 3.5 K; 
solid lines connect bands which are produced by similar l attice 
components; dashed lines indicate fits obtained in Section 6.1.1 



Some bands are missing from the forsterite and enstatite spectra 
because the samples were either too thick, or too thin for them 
to be resolved. 



5.4 Wavelength Shift and Splitting Patterns 

In Figure ^ the shifts in the peak wavelength of the bands 
between 295 K and 3.5 K, A(295 K)-A(3.5 K), are plotted as 
a function of the wavelength of the respective 295 K peaks. 
In general, bands shift to shorter wavelengths at 3.5 K and 
the magnitude of the band shift increases with increasing 
wavelength, from ~ — 0.2 /xm at 25 /xm to ~ 0.9 /^m at 
80 fim.. However, the shift patterns are quite complicated 
because different peaks have different origins in the lattice. 



The situation is complicated in some of these spectra since 
we have not been able to resolve all the peaks in the thick 
forsterite and thin enstatite. Measurements of the spectra 
of other film thicknesses are required to determine the low- 
temperature behaviour of all the infrared peaks. 

The shift patterns and band splittings are similar within 
a solid solution series if the lattice structure is the same. 
For example, the ~ 50-/im and ~ 70-/im bands exhibited by 
forsterite and olivine shift by similar amounts, -0.5 pm and 
-0.9 /xm, respectively, when the minerals are cooled from 
295 K to 3.5 K. The increase in the energies of these Si04~ 
and combined Si04~ and Mg'^"'' translations at 3.5 K in- 
dicates a contraction of the lattice as the temperature is 
reduced; sharpening of the bands is due to a reduction in 
phonon scattering at 3.5 K. 

The ortho- and clinopyroxenes exhibit different shift 
and splitting patterns. The band shifts in Bramble en- 
statite and in enstatite increase fairly uniformly with wave- 
length. However shifts between neighbouring bands in diop- 
side clinopyroxene vary by up to 0.3 /xm and between 27 and 
41 /im the pattern is almost saw-tooth. The complex shift 
pattern of diopside is probably caused by the difference in 
the size and mass of the Ca'^"'' and Mg^^ cations which re- 
spectively occupy the M2 and Ml sites. Band splitting in the 
three 18-22. 5-/im bands is common to both orthoenstatites- 
the Si-O-Si shifts in both minerals are ~ 0.0 and -0.4 
and the Mg^""" shifts of the ^ 22.5-^m peak are 0.0 and 
-0.5 /xm. The splittings are caused by a reduction in the 
widths of the peaks, and by different orientations (polariza- 
tions) of the bonds within the lattice because the different 
crystallographic axes contract with temperature at different 
rates. 



6 DISCUSSION 

6.1 Astrophysical Implications 

These determinations of the relationship between tempera- 
ture and band shifts and band splitting could help in reach- 
ing a better understanding of the mineralogy of astronom- 
ical regions showing optically thin silicate bands in the far 
infrared. Such temperature-dependent laboratory measure- 
ments could greatly benefit the elucidation of the physical 
properties of dusty environments. 



6.1.1 Limits on the effect of temperature on Xpk 

The far-infrared bands observed in many astronomical en- 
vironments are wavelength-shifted with respect to bands 
of similar shape observed in room-temperature laboratory 
spectra (Molster 2000). Wavelength shifts can be intro- 
duced by temperature differences, but they can can also 
be introduced by other factors such as chemical and iso- 
topic composition and degree of crystallinity. To place lim- 
its on the contribution of temperature to the wavelength 
shift of a peak in astronomical data, fits to the ratio, 
R = (A(295K) - A(3.5K))/A(295K), for peaks which did not 
split at 3.5K were obtained over the wavelength ranges spec- 
ified in Table ^. Linear fits to R were obtained using the IDL 
LINFIT routine; the fit to the upper limit of the diopside 
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Table 3. Fits to Rx = A{295K) - A(3.5K) using polynomial fits 
of the form, R = a + bx + cx^, where x=Apfe(295 K), for the 
wavelength shift data plotted in Figure g. 



Mineral 


Range 


a 


b 


c 




/im 


xlO-3 


xlO-3 


xlO-3 


Forsterite 


31-70 


-9.66 


0.347 




Olivine 


24-74 


7.47 


0.0396 




Enstatite 


24-44 


-14.0 


0.686 




Bramble Enstatite 


23-52 


-3.13 


0.356 




Diopside Lrl 


19-67 


-5.44 


0.262 






25-45 


62.4 


-2.55 


0.0294 



" Lower boundary of wavelength shifts: fit to 295 K peaks at 
19.38, 23.17, 27.19, 32.30 ^im and all peaks where A > 40 ^m. 
' Upper boundary of wavelength shifts: fit to 295 K peaks at 
25.22, 34.29, 44.93 

shifts was obtained with the IDL POLYFITW routine. Poly- 
nomial fits to the wavelength shifts of singlet bands between 
295K and 3.5K are plotted in Figure ^ and the coefficients 
of these fits are are listed in Table ^ 

6.1.2 Determination of grain temperatures 

Our results indicate that suitable high resolution spectra ob- 
tained at a range of cryogenic temperatures can be used to 
determine the characteristic temperatures of specific grain 
populations in astronomical spectra. An emission feature 
which peaks near 69 /xm in circumstellar spectra has been 
identified with forsterite (Malfait et al. 1998). However, the 
astronomical feature is shifted relative to room temperature 
laboratory data, which give Xpk = 69.74 /xm (Figure |^ and 
Table In Figure ^b) we compare the shape of the 69- /im 
band of our laboratory forsterite spectra obtained at three 
different temperatures. The spectral resolution of the ISO 
LWS grating spectrometer at this wavelength was 0.3/xm, 
compared to the laboratory feature's FWHM of 1.3 fj,m at 
295 K and 0.4 fim at 3.5 K, while the LWS wavelength preci- 
sion was 0.03/im, compared to the shift of 0.9/im in the peak 
wavelength between 295 K and 3.5 K. Therefore the deriva- 
tion of characteristic grain temperatures from ISO spectra is 
quite feasible (see also Molster 2000) and will be reported on 
in a subsequent paper (Bowey et al., in preparation). Band 
splitting and band shifts at other wavelengths, for example 
the splitting in the 36.5/im forsterite band shown in Figure 
^a), could also be used to determine grain temperatures if 
the spectral resolution and signal-to-noise of the astronom- 
ical data were sufficiently high. The relation between the 
derived characteristic grain temperatures at different wave- 
lengths may aid the determination of the spatial distribution 
of the dust. 

6.1.3 Effect on inferred mineralogical compositions and 
dust masses 

The absorption coefficients of the forsterite bands are en- 
hanced at 3.5 K compared to room temperature values. 
In contrast, there is little change in the absorption coeffi- 
cients of most of the pyroxene features. Therefore, the true 
mass ratio of forsterite to enstatite could be smaller than 
estimated using room-temperature laboratory data. More 
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Figure 7. Comparison between 3.5-K (solid), 77 ± 5-K{dashed) 
and 295-K (dotted) forsterite spectra). 77 K spectra were obtained 
by placing the forsterite sample in a cold finger cooled by liquid 
nitrogen within the output beam of the spectrometer. 

mass is required to reproduce the strengths of many fea- 
tures in ISO spectra using room temperature optical con- 
stants, than would be required at the likely temperatures of 
the astronomical grains. Measurements of the temperature- 
dependence of the mass absorption coefficients are required 
in order to obtain the true relative mineralogical composi- 
tion and dust masses. 



7 SUMMARY 

We have obtained 0.25 cm~^ resolution 15-85 fj,m spec- 
tra of powdered minerals embedded in petroleum jelly 
at 295 ± 2 K and 3.5 ± 0.5 K. The samples were two 
olivines, forsterite (Mg2Si04) and olivine var. peridot 
Mgi,77Feo.2oAlo.oiNio.oiSi04, two orthoenstatites, enstatite 
and Bramble enstatite, and a clinopyroxene (MgCa diop- 
side). 

In general: 

(i) the room temperature spectra look like lower- 
resolution versions of their 3.5 K counterparts. In both of 
the olivines the degenerate ~ 21 fim Mg^"*" translation and 
Si04~ rotation splits into two bands when cooled to 3.5 K. 
The ~ 34-^m peak in forsterite is saturated, however this 
Mg^+ translation band is split in the 3.5 K olivine spectrum 
and the 36.45-^m 295 K Mg^^ translation peak in forsterite 
is also split at 3.5 K. In the orthoenstatites, band splitting 
is seen only in the three bands shortward of 23 fim. In the 
MgCa diopside clinopyroxene only the 20. 91- fim degenerate 
Si-O-Si deformation is split. Different sample thicknesses 
must be measured to verify that no splitting is seen in other 
bands. 

(ii) the wavelength of a peak decreases with temperature 
and the magnitude of the shift increases with the wavelength 
of the band, from ~ — 0.2/im at 25 fj,m to 0.9 pm at 
80 j-im. However, the shift patterns are quite complicated 
and depend on lattice structure as well as stoichiometry 
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since they differ between ortho- and clino- pyroxenes. For 
example, the band shifts in Bramble enstatite and enstatite 
increase with wavelength. However, in the diopside clinopy- 
roxene the shifts between neighbouring bands vary by up to 
0.3 i-im - between 27 and 41 /xm the pattern is almost saw- 
tooth. The situation is complicated in these spectra since 
we have not been able to resolve all the peaks at all tem- 
peratures. For the olivines and enstatites the wavelength 
shifts of bands which are singlets at 295K and 3.5K can be 
approximated by polynomials of the form ax + bx^ where 
X = Apfe(295 K) and coefficients a and b differ between min- 
erals. Lower limits to the wavelength shifts for diopside can 
be obtained using polynomials of this form, however a third 
order term is required to obtain the upper estimates for 
wavelength shifts in the 20-41/im region. 

(iii) bands sharpen as the temperature is reduced. At 
3.5 K the widths of isolated single bands in the olivine, 
enstatites and diopside are 90% of their 295 K-widths. 
However, in forsterite, the 3.5-K, 31, 49 and 69-pm band- 
widths are, respectively, 90%, 45% and 31% of their 295 K 
widths. 

(iv) the relative absorbances of neighbouring peaks 
change and bands sharpen as the temperature is reduced; 
this effect is most noticeable in the forsterite spectra. Peaks 
at 53.91 fim in Bramble enstatite and at 57.16 fj,m in diop- 
side appear only in the 3.5 K spectra. 

Our results show that comparison between the band 
shifts, bandwidths and band splitting in laboratory data 
obtained at a range of temperatures with astronomical far- 
infrared spectra can allow a determination of the character- 
istic temperatures of crystalline grains responsible for the 
optically thin emission features. Further low-temperature 
spectra at a range of temperatures are required to ascer- 
tain the precise temperature dependence of the shift, since 
the relationship between T and Xpt is unlikely to be linear 
between 295K and 3.5K and it is possible that the bands 
stop shifting at some temperature above 3.5K. Ultimately 
it may be possible to determine the temperature of distinct 
grain populations by comparison of the observed peak wave- 
lengths and FWHMs with laboratory spectra. Since the ab- 
sorption coefficients of bands in laboratory spectra also vary 
with temperature, new measurements may allow better esti- 
mations of grain masses and silicate mineralogy. Additional 
high-resolution laboratory spectra of a larger variety sili- 
cates, at a larger number of temperatures, are desirable for 
comparison with ISO data and with future measurements 
obtained by spectrometers on board SIRTF, FIRST and the 
NGST. These laboratory measurements are planned. 
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